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Abstract 

We discuss neutral Higgs boson production through gluon fusion and bottom-quark 
annihilation in the CP-conserving Z 3 -invariant Next-to-Minimal Supersymmetric Stan¬ 
dard Model (NMSSM) at proton colliders. For gluon fusion we adapt well-known asymp¬ 
totic expansions in supersymmetric particles for the inclusion of next-to-leading order 
contributions of squarks and gluinos from the Minimal Supersymmetric Standard Model 
(MSSM) and include electroweak corrections involving light quarks. Together with the 
resummation of higher-order sbottom contributions in the bottom-quark Yukawa cou¬ 
pling for both production processes we thus present accurate cross section predictions 
implemented in a new release of the code SusHi. We elaborate on the new features of 
an additional SU(2 )l singlet in the production of CP-even and -odd Higgs bosons with 
respect to the MSSM and include a short discussion of theoretical uncertainties. 
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1 Introduction 


After the discovery of a scalar boson at the Large Hadron Collider (LHC) [T|[2] in 2012 an essen¬ 
tial task of particle physicists is to reveal the nature of the Higgs-like state and thus the nature 
of electroweak symmetry breaking. Apart from deviations from the Standard Model (SM) 
prediction of the properties of the found Higgs-like state, further work includes the search 
for additional less and/or more massive scalar bosons, which can nicely be accommodated in 
supersymmetric models. The Next-to-Minimal Supersymmetric Standard Model (NMSSM) 
extends the Minimal Supersymmetric Standard Model (MSSM) by an SU(2)/, singlet and 
allows the dynamical generation of the //-term through electroweak symmetry breaking W- 
The latter singlet-doublet mixing term in the superpotential lifts the MSSM tree-level upper 
bound of the Higgs mass given by the Z-boson mass. Thus, the NMSSM can easily accommo¬ 
date the SM-like Higgs boson with a mass close to 125 GeV. Whereas for the calculation of the 
NMSSM Higgs spectrum and branching ratios various spectrum generators are available and 
include higher orders in perturbation theory (see Section [3]), the calculation of neutral Higgs 
production cross sections did not exceed leading order (LO) in quantum chromodynamics 
involving squarks and gluinos (SQCD) [5] and did not include electroweak corrections - apart 
from private implementations in e.g. HIGLU [6]. 

It is therefore timely to present the missing ingredients and a code for the calculation of 
accurate neutral Higgs production cross sections in the NMSSM, where the five neutral Higgs 
bosons are predominantly generated through gluon fusion and bottom-quark annihilation at 
a proton collider. For this purpose we extend the code SusHi [7]. For the time being we 
restrict our implementation to the real NMSSM without additional CP violation, such that 
CP-even H\,H 2 and H$ and CP-odd Higgs bosons A\ and A 2 can be distinguished in the 
Higgs sector. Most recent efforts related to Higgs physics at the LHC are summarized in 
the reports of the LHC Higgs cross section working group [MO]. The SM Higgs is mainly 
produced through gluon fusion, where the Higgs-gluon coupling is mediated through virtual 
top- and bottom-quarks m Higher order QCD corrections at next-to-leading order (NLO) 
are of large importance P2HH]. In the effective theory of a heavy top-quark the inclusive 
cross section is known to next-to-next-to-leading order (NNLO) in QCD [151117] . in addition 
finite top-quark mass effects at NNLO were calculated 1 181122] . Beyond NNLO QCD effects 
are accessible through resummation I23H28] and electroweak corrections are known (291131] . 
Meanwhile next-to-NNLO (NNNLO) QCD contributions were estimated in the so-called thresh¬ 
old expansion (321435] . but they are not further considered in this publication. 

The SM results for Higgs production through gluon fusion can be adjusted to the MSSM 
and the NMSSM through a proper reweighting of the Higgs couplings to quarks. However, the 
gluon fusion process can also be mediated through their superpartners, the squarks. With 
respect to the MSSM the only generically new ingredient, which goes beyond the projection 
of the physical Higgs bosons onto the neutral components of the two Higgs doublets, are 
couplings of the NMSSM singlet to squarks, since no couplings of the singlet to quarks or 
gauge bosons are present in the tree-level Lagrangian. It is therefore of importance to include 
squark contributions to gluon fusion at the highest order possible, even though they decrease 
in size with increasing squark masses. For the pseudoscalars A, squark contributions to gluon 
fusion are only induced at NLO, which motivates to go beyond just LO squark contributions 
for all Higgs bosons. For this purpose we adapt the works of Refs. [35H58] for the MSSM to 
present NLO SQCD contributions for the NMSSM, which are based on an expansion in terms 
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of heavy supersymmetric particles taking into account terms up to 0(m|/M 2 ), O(rrif /M 2 ), 
0(m 2 /M 2 ) and 0(m 2 z / M 2 ), where denotes the Higgs mass and M a generic SUSY mass. 
In contrast to the MSSM we are at present only working in this expansion of inverse SUSY 
masses and do not include an expansion in the so-called VHML, the vanishing Higgs mass 
limit {rriff, —>• 0) for the SQCD contributions, as implemented in evalcsusy j39ff41] or discussed 
in Ref. [32] • In the latter limit higher-order stop-induced contributions up to NNLO level are 
known [431445] and were partially included in previous discussions of precise MSSM neutral 
Higgs production cross sections [46] . Although for a pure CP-odd singlet component NNLO 
stop-induced contributions are the first non-vanishing contributions to the gluon fusion cross 
section, we leave an inclusion of these to future work. For completeness, we add that in the 
MSSM a numerical evaluation of NLO squai'k/quark/gluino contributions was also reported in 
Ref. [471138] , whereas Refs. [39U5I] presented analytic results for the pure squark induced NLO 
contributions. Electro-weak contributions to the gluon fusion production process mediated 
through light quarks 30.31] can be adjusted from the SM to the MSSM 7] and similarly to 
the NMSSM and are known to capture the dominant fraction of electroweak contributions for 
a light SM-like Higgs with a mass below the top-quark mass, whereas they are generically 
small for larger Higgs masses. 

For large values of tan f3 , the ratio of the vacuum expectation values of the neutral com¬ 
ponents of the two Higgs doublets, the bottom-quark Yukawa coupling is enhanced, such that 
the bottom sector gets more important for gluon fusion, and the associated production with 
a pair of bottom-quarks pp —>• bb(j) is significantly enhanced. SusHi includes bottom-quark 
annihilation bb —» (f>, which in case of non-tagged final state 6-quarks is a good theoretical 
approach, since it resums logarithms through the 6-parton distribution functions. The latter 
process is known as five-flavor scheme (5FS) up to NNLO QCD [52, 53] and can easily be 
reweighted from the SM to the MSSM/NMSSM by effective couplings [531155] , In the NMSSM 
the singlet does not couple to quarks at LO, however taking into account the singlet induced 
component into the resummation of higher-order sbottom effects is mandatory, since also the 
singlet to sbottom couplings are enhanced by tan (3. 

The new release of SusHi thus provides gluon fusion cross sections at NLO QCD taking 
into account the third generation quarks and their superpartners, the squarks, for all the 
five neutral Higgs bosons of the NMSSM. The squark and squark/quark/gluino contributions 
are implemented in asymptotic expansions of heavy SUSY masses. Electro-weak corrections 
induced by light quarks through the couplings of the Higgs bosons to Z and W ± bosons can be 
added consistently like in the MSSM. Similarly, the NNLO top-quark induced contributions are 
included. In addition, sbottom contributions can be resummed into an effective bottom-quark 
Yukawa coupling, also taking into account the additional singlet to sbottom couplings. The 
latter also applies to the calculated bottom-quark annihilation cross section at NNLO QCD. 
All features SusHi provides for the MSSM are available for the NMSSM as well, in particular 
distributions with respect to the (pseudo)rapidity and transverse momentum of the Higgs 
boson under consideration can be obtained. Left for future work is a link to MoRe-SusHi [56] 
to allow for the calculation of momentum resummed transverse momentum distributions. 

We proceed as follows: We start with a discussion of the theory background in Section [21 
where we elaborate on the NMSSM Higgs sector and the calculation of the gluon fusion cross 
section. Then we present the NLO virtual amplitude for gluon fusion as well as the calculation 
of bottom-quark annihilation including the resummation of sbottom-induced contributions 
to the bottom-quark Yukawa coupling in the NMSSM. Subsequently we comment on the 
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implementation in SusHi in Section [3l before we investigate the phenomenological features of 
the singlet-like Higgs boson in the CP-even and CP-odd sector with regard to Higgs production 
in Section [4) We also include a short discussion of theoretical uncertainties. Finally, we 
conclude and present the Higgs-squark-squark couplings in Appendix El 

2 Theory background 

In this section we discuss the Higgs sector of the CP-conserving Z 3 -invariant NMSSM, before 
we proceed to the resummation of tan/3 enhanced sbottom contributions in the bottom- 
quark Yukawa coupling. Subsequently we move to the discussion of the Higgs production 
cross section in gluon fusion, where we present the adapted formulas for the NLO SQCD 
virtual amplitude, and finally comment on the consequences of the additional singlet to 
bottom-quark annihilation. 

2.1 The Higgs sector of the CP-conserving NMSSM 

Our notation of the Higgs sector of the CP-conserving Z 3 -invariant NMSSM closely follows 
Ref. [57]. For NMSSM reviews we refer to Refs. [31IU. The superpotential can be written in 
the form 

ITnmssm = Wmssm - e ab \SH a d H h u + -kS 3 , (1) 

where Wmssm equals the superpotential of the MSSM without //-term. S denotes the ad¬ 
ditional SU(2 )l singlet superfield compared to the MSSM with the two SU(2)^ doublet su¬ 
perfields H d and H u . e ab contracts the SU(2)^ doublet components. Since the singlet S is 
a neutral field, it induces one additional CP-even and one additional CP-odd neutral Higgs 
boson as well as one additional neutralino compared to the MSSM. The soft-breaking terms 
include the scalar components H d , H u and S of the superfields and are given by 

£soft = £ so ft,MSSM + (eab^A\SH%H b — -kA k S 3 + h.c.) — m 2 s \S\ 2 . (2) 

The soft-breaking mass m s can be derived from the minimization conditions of the tadpole 
equations (in addition to ni 2 Hd and m 2 Hu like in the MSSM), whereas A\ and A K are usu¬ 
ally considered input parameters. A\ can be alternatively replaced by the charged Higgs 
mass riiH± as input parameter. The neutral components of the Higgs fields are decomposed 
according to 

H * = fy Vd + H d+iH I d ), H^ = ^=(v u + H^ + iHi), S = ^=(v s + S R + iS I ), (3) 

where v d , v u and v s denote the vacuum expectation values (VEVs) and the fields with indices 
R and / are the CP-even and CP-odd fluctuations around them. An effective term is 
generated through the VEV of the singlet 


which will be further used within this article. We do not present the explicit form of the 
mass matrices here, but refer to Ref. m- We define the CP-even gauge eigenstate basis 
H r = (H R , H R , S R ) and the CP-odd one H 1 = (H^, H^, S 1 ). Whereas in the former case 
the mass eigenstates Hi with i € {1, 2, 3} are obtained through one rotation 


3 = 1 


R 


(5) 


we perform a prerotation in the CP-odd sector to obtain the MSSM pseudoscalar A and the 

Goldstone G in the form H' 1 = (G, A, S 1 ). The prerotation 1Z G is given by the ratio of 

vacuum expectation values tan/3 = v u /vd, such that the mass eigenstates Ai with i € {1,2} 
are obtained by 

33 /eg -sp 0\ 

A l = J 2 ^f+i,j H j I =J 2 7Z ^ 3 n fkH I k with K G =\sp cp 0 , (6) 

j= 1 j,k= 1 \0 0 1/ 

where TZ P is a (3 x 3)-matrix, which however only consists of a (2 x 2)-mixing block, whereas 
7Z r = 7Z r = 0 for i / 1 and TZfl = 1. For Higgs production the Goldstone boson does 
not need to be considered. In the following we make use of the notation “singlet-like Higgs 
boson”, which refers to the CP-even/odd Higgs boson with the dominant fraction of the 
singlet component S in gauge eigenstates. For this purpose we define the singlet character 
|Ff 3 | 2 for Hi and \H p +l ^\ 2 for A % . In our discussion of cross sections we denote the Higgs 
boson by the letter i/>, which can be replaced by any of the physical Higgs bosons Hi or Ai. 
For picking viable scenarios for phenomenological studies we refer to Ref. [58] for a recipe to 
obtain positive eigenvalues for the singlet-like CP-even and -odd Higgs by varying A K between 
a minimal and a maximal value. 

Whereas the singlet component S does not couple to quarks, F-terms induce a coupling of 
the singlet-like Higgs to squarks, which is of relevance for Higgs production. We present the 
Higgs-squark-squark couplings to the third generation of squarks in Appendix [A] We point 
out that the singlet component mixes with the Higgs doublets proportional to A and also 
the couplings of the singlet component to squarks are proportional to A. It is thus possible 
to mostly decouple the singlet component by lowering the value of the parameter A. The 
couplings to quarks can be easily translated from the MSSM by the correct projection on the 
neutral doublet components H R , H R and the pseudoscalar A and yield relative to the SM 


Hi rr> S 1 


ll cos /3 ’ 


9? = K 
9d 1 = K + l,2 


tan/3, 


9u l = K? 2 


9u 1 =K 


sin /3 


■i+1,2 


tan/3 


(7) 


with i € {1,2,3} in the CP-even and i € {1,2} in the CP-odd Higgs sector. The relative 
strength g ^ enters the Yukawa couplings in the form Yj = \/2m.fg^/v with the vacuum 
expectation value v 2 = + v 2 . 


2.2 Resummation of higher-order sbottom contributions 

It is well-known in the MSSM that tan /3 enhanced sbottom corrections to the bottom-quark 
Yukawa coupling can be treated in an effective Lagrangian approach [591464] to be resummed. 
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For the case of the NMSSM just taking into account SQCD corrections the effective Lagrangian 
can be written in the form |65] 


-Ceff = —YbbR 


2 a. 


H°,+ 


AA b 
[i tan (3 


S*H\ 


o* 


b L 


with 


Ab = -—ms^tan/3/(m~ , rnf ,m~) and I(a,b,c ) 
0 7r t>i b 2 u 


ab log l + bc\og- c + ca log § 
(a — b) (b — c)(a — c) 


( 8 ) 

• (9) 


Ref. [65] additionally presents the inclusion of SUSY electroweak corrections proportional to 
the soft-breaking parameter At- The inclusion of electroweak corrections into does not 
harm our subsequent discussion of SQCD corrections and can thus always be included in the 
bottom-quark Yukawa coupling entering gluon fusion and bottom-quark annihilation. Apart 
from a coupling of the bottom-quarks to the gauge eigenstate H® also an effective coupling to 
the singlet S can be induced at loop level, the latter being proportional to \v u instead of p. 
The sbottom corrections can be absorbed into effective Yukawa couplings, which read [65] 


- Hi 

9 b = 


Hi 

9 b 


1 + A b 


1 + A b 

for the three CP-even Higgs field Hi and 


2 


+ 


Hi 


i3 vcos /3 


tan /3 TZf 3 v 


~Ai 

9 b = 


Ai 

9b 


1 + Ab 


1 + Ab — 


H'i+l,3 V 


tan 2 (3 2 v s tan (3 


( 10 ) 


( 11 ) 


for the two CP-odd Higgs fields Ai. 


2.3 Gluon fusion cross section 


After our discussion of the CP-conserving NMSSM Higgs sector and the resummation of 
sbottom contributions in the bottom-quark Yukawa coupling, we present the gluon fusion 
production cross section for a Higgs boson 0, which can be written in the form [7] 


a(pp -A 4> + X) = a 4 Q 


l + C* 


a. 


7 r J 


-b A<jg g + A at + A at 


dTH, 


gq 


qg 


( 12 ) 


with = rn^/s and the hadronic centre-of-mass energy s. The factor Uq includes the LO 
partonic cross section. encodes NLO terms of singular nature in the limit s -A m^ with 
the partonic centre-of-mass energy s. The gluon-gluon luminosity is given by the integral 

dC" f 1 dx 

— = l -0( X )s(r/ X ) ■ (13) 


The contributions Aaf g , Aaf q , and A<t)T are the regular terms in the limit s -A in 
the partonic cross section and arise from gg, gq and qq scattering, respectively. The LO 
contribution ay is obtained by the formulas presented in Ref. |7| using the NMSSM couplings 
of the Higgs bosons to quarks and squarks. Similarly the contributions A a£ y are obtained 
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from the MSSM by a proper replacement of the involved Higgs boson to quark and squark 
couplings. The factor C& can be decomposed in the form 

+ 7T 2 + /3o log fe^ , (14) 

with (3q = 11/2 — rif/3 and ny = 5 as well as the factorization and renormalization scales, 
/i F and \i R respectively. is the LO (one-loop) virtual amplitude in the limit of large 

stop and sbottom masses. is the NLO (two-loop) virtual amplitude and equals the form 
factors ; for the CP-even Higgs bosons and A) for the CP-odd Higgs bosons, which are 
presented in the following Section 12.41 to account for NLO virtual contributions from quarks 
and squarks in an appropriate way. In the MSSM limit they correspond to the form factors 
of Refs. [364438] except from a constant factor of —3/4. Contrary to the case of the MSSM 
we do not employ evalcsusy |401l4lj to obtain the NLO amplitude in the limit of heavy 
top-quark and stop masses for the light Higgs, but use the expanded form factors presented 
in the following sections instead. Accordingly our implementation does not (yet) include 
approximate NNLO stop contributions as presented in Ref. US] for the MSSM. The NNLO 
top-quark contributions in the heavy top-quark effective theory making use of Refs. [151166] 
are included according to Eq. (29) of Ref. [7j. 

Lastly we comment on the inclusion of the electroweak corrections to the gluon fusion 
production cross section. Similarly to the MSSM the full SM NLO electroweak (EW) correc¬ 
tions (23] can be added to the top-quark induced result only, assuming complete factorization 
of EW and QCD effects EZj. We recommend the latter procedure only for a SM-like Higgs 
boson. Contrary the inclusion of electroweak corrections due to light quarks ISME], where 
the Higgs boson couples to either the Z or boson, can be adjusted to the MSSM and 
accordingly the NMSSM in an appropriate way [68] • For this purpose the generalized couplings 

ffv ‘ = ftfi cos /3 + sin fj (15) 

of the i-th CP-even NMSSM Higgs boson to the heavy gauge boson V € {W^jZj need to 
be inserted in the formulas of Ref. [7J. The missing projection 7 Zf 3 on the singlet component 
reflects the fact that the singlet does not couple to gauge bosons. The CP-odd Higgs bosons 
do not couple to gauge bosons either, such that electroweak corrections due to light quarks 
are absent. 

2.4 NLO virtual amplitude for gluon fusion 

Regarding the implementation of two loop contributions to gluon fusion, we closely follow 
Refs. [364438] for the MSSM, which can be translated to the NMSSM. Their calculation at NLO 
is based on an asymptotic expansion in the masses of the supersymmetric particles. We can 
project the form factors onto the ones in gauge eigenstates according to 

uf l) = ~(nf A uf {2l) + nf 2 u^ 2l) + nis R ^ 2l) ) (16) 

A (2 ° = + < +1 , 3 S J ’ (a) ) • (17) 

The individual contributions in gauge eigenstates are presented in Section 12.4.11 for the CP- 
even and in Section 12.4.21 for the CP-odd Higgs bosons. We included the constant factor 
between Refs. [364138] and our work in the above equations. 


C* = 2Re 


( j >( 2 0 

<j>dO 
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2.4.1 CP-even Higgs bosons 


In this subsection we present the form factors for the CP-even Higgs bosons in gauge eigen- 
state^] 


R,(2l) _ 


K 

njR,(2l) 

' L u 

$R,( 20 = 


1 


-m t gs 2 Q t F? 1 + m z s 2/3 D't 


sin/3 

1 

cos/3 . 

1 


21 


-rn b l-is 20b F t f - m z s 2 pD b 


+ 


+ 


1 


m b A b s 2 o b Fg l + 2 m b Gf + 2 rri z z CpD b 


2r>21 


,2 2 r>2l 


sin /3 


1 


--^ m A VdS 2e t F t 


21 


+ 


cos /3 


cos /3 

1 

sin/3 . 

-±=m b \v u s 26b Ft 


m t A t s 2 e t Ff + 2 mfG'f - 2m z s Z pD z 


2ri2l 


■2 2 r\2l 


21 


(18) 


which includes the effective g parameter defined in Eq. Q. All functions in 'Hu 

and S R, ( 2 0 can be directly taken over from Refs. [36(138]. keeping in mind the different 
convention in the sign of the /i parameter. For on-shell (OS) parameters (see Refs. [361138] 1 
and thus for our implementation the contribution Ff l * is shifted according to Section 3.3 of 
Ref. [38J and i 7 b l according to Ref. [36]. The shift also applies to F b l and F^ 1 entering the 
singlet contribution S R , since the differences in the prefactors being g, Xv^ or Xv u are not 
renormalized when taking into account SQCD contributions and therefore do not contribute 
to the described OS shifts. 

It remains to discuss the inclusion of resummed sbottom contributions into the bottom- 
quark Yukawa coupling within the virtual corrections to gluon fusion, where care has to 
be taken to avoid a double-counting of NLO SQCD contributions. The naive resummation 
9b = 9b/( 1 + A;,) is incorporated in the same way as in case of the MSSM |36y36]. The 
resummation as presented in Section 12.21 instead needs the subtraction of the tan /3 enhanced 
contributions to G'f multiplied with the corresponding coupling correction, in detail for the 
three CP-even Higgs bosons Hi 


2 m 2 b Gl l -> 2 m 2 b Gl l - ^yGi /2 (r b )fi tan£ m^, mj)) K[ 

with the factor K | being 
1 


K' = 


1 + A b 


1 - 


K?2 

TZf x tan /3 


+ 


v cos /3 

Ff x v s 


*€{1,2,3} 


(19) 


( 20 ) 


All occurrences of the bottom-quark Yukawa coupling in the two-loop amplitude are multi¬ 
plied with the factor Ki = g Ri / g Ri using Eq. ([ID]), such that the shift reported in Eq. (fTDl) 
avoids double-counting of the purely SQCD induced contributions at the two-loop level. Em¬ 
ploying the expansion in heavy SUSY masses the NLO virtual contributions to neutral CP-even 
Higgs production in the NMSSM are now fully presented. 


2.4.2 CP-odd Higgs bosons 

We now turn to the case of the two CP-odd Higgs bosons, where we present the form factor 
in the basis H’ 1 after a prerotation from gauge eigenstates 1 . At LO only diagrams involving 

1 In the CP-even sector we adapt the MSSM results of Refs. I36U38) to the NMSSM by isolating the terms 

proportional to the Jl^/Jl^-squark-squark couplings and replacing them by the S fl -squark-squark couplings 

for the form factor S. Similarly we proceed in the CP-odd sector starting from the form factors of Ref. m 

taking into account the prerotation of the CP-odd Higgs mixing matrix. 
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quarks coupling to the pseudoscalar A exist, such that the form factor at LO only consists 
of the part whereas S equals zero. At NLO however couplings of A and S 1 to 

squarks qiqj for i / j induce contributions to Higgs production. The two-loop form factor 
presented in Ref. m for the MSSM can therefore be translated to 


n 


1 ,( 21 ) _ 


cot (3(lCf g + /C 2 ~~) + tan /3(lCf g + /C 2 ~ ) 


-21 
"tig 


21 


-21 
" bbg> 


S / ’ (2Z) = [cot [tK S f + tan /3 TCf^ 

' ttg ' bbg 


( 21 ) 


Whereas the individual contributions to can be taken from Ref. m, we present the 

contributions to S 1 ^ 21 * separately: 


JC 


S,2l_C Flr u f ' m 

ttg 


~>C U (rt) 


g m *72 Xv 


2 v L ' nit m~ — m 2 


*2 


x\ 


4 


1 — x\ 


In x\ — -——- In x\ 


m t „,t . 
— y—S26» i 'V l + 


rrtn 


-S,2i _C F n, ,m,g m b ^j 


2m t ^2 Xv 


i 


1-4 


n 2 -^-ic ll (T t ) 


rn % 


rria(m2 — m? ) 2 7 m~ — m 2 

y V f-, C2 7 


Cl 


-7^4 


^2 


iT~ cV/ \^b)~ o o 

bb 9 2 m b mi — m~ 


1 — 4 


In 4 — 


1-4 


In 4 ) — ——' 


s 2 e b 7Z!\ 


( 22 ) 

(23) 


bl b 2 

with x\ = m 2 /m~ and 4 = m~ /m~. The functions 1C 11 and 1Z2 can be found in Ref. m- 
The functions 7Z'\ and 77 /\ are given by: 


1Z’\ = 


tz!\ = 


c F 


7=2 Xv 


(x\ — 4) 2 m g 
— (x\ •H' x\) 

C F 


x +2 


4 2 (1 — 24 ) x\(xY — 2x1, + x^Xg) lnx} 


[ 2(1 — 4)(i — 4) 


+ 


2(1 - x\) 2 


(24) 


75 Au 


(x\ — Xg) 2 m.g 


x[ 2 (l — 2^2) x\(xi — 2 x\ + X4X2) lnx| 


2(! — 4)(! — 4) 


+ 


2(1 - x { ) 2 


In the bottom sector the relevant function yields: 


TZf = 


C F 


7^ Xv 


Xi(l — 2xg) x±(xi — 2 x ° 2 + X4X2) lnx* 


(4 - 4) 2 m 9 L2(i-4)(!-4) 


+ 


2(1 - x’5) 2 


- (; x\ 4) (25) 


-(x44) (26) 


The shifts of individual contributions in case of OS parameters can be taken over from the 
MSSM case. The inclusion of resummed sbottom contributions to the bottom-quark Yukawa 
coupling needs the following shift in the two-loop form factor for the CP-odd Higgs bosons At 

C F 




^blg ~ -y7C lz (r b )/r tan/3 (-rrigliml^ml^m 2 ^ K[ 


using 


K[ = 

1 1 + A b 


1 + 


+ 


n i+ 1,3 V 


tan 2 (3 7 Z^ +1 2 v s tan /3 


*€{ 1 , 2 } 


(27) 


(28) 


Again we point out that our sign convention with respect to g is opposite to Ref. m and all 
occurrences of the bottom-quark Yukawa coupling in the two loop amplitude are multiplied 
with the factor I<i = using Eq. (fill) . 


































2.5 Bottom-quark annihilation cross section in the 5FS 

The generalization of the calculation of bottom-quark annihilation cross sections in the five- 
flavor scheme (5FS) from the MSSM to the NMSSM case is straightforward by using the 
appropriate couplings of Higgs bosons to bottom-quarks. For this purpose the resummation 
of sbottom contributions as described in Section 12.21 is taken into account. For the specific 
case of the singlet-like Higgs boson we point out that in case the coupling to the bottom-quark 
vanishes (due to cancellations in the mixing with the Higgs doublets) a priori the coupling 
to sbottom squarks can still be present. This is not taken into account by the resummation 
procedure. 


3 Implementation in SusHi 


In the current implementation of neutral Higgs production in the real NMSSM within the 
code SusHi the Higgs mixing matrices as well as the Higgs masses have to be provided as 
input in SUSY Les Houches Accord (SLHA) form [691170] and can be obtained by spectrum 
generators for the NMSSM. Common codes are NMSSMTools [ZHEH, NMSSMCALC [57U6M75B T7] . 
SOFTSUSY [781179] . SPheno+Sarah [801182] and FlexibleSUSY+Sarah [801153] , 

Special attention needs to be paid to the renormalization of the stop and sbottom sector, 
which in the ideal form should be identical in the calculation of Higgs masses and mixing and 
the calculation of Higgs production cross sections. For the time being, SusHi either relies on 
the internal calculation of on-shell stop and sbottom sectors as described in the manual [7] 
or on the specification of the on-shell masses m and rriq 2 and mixing angles Oq in the input 
file. For both cases input files can be found in the folder /example within the SusHi tarball. 
The user is asked to check the meaning of output parameters of spectrum generators, i.e. the 
chosen renormalization scheme. If the user specifies the on-shell squark masses and mixing 
angles together with the on-shell soft-breaking parameters At and Af } by hand, she/he should 
make sure that in the stop sector At as well as the on-shell top-quark mass mt, the on-shell 
stop masses mi and m^ 0 and the mixing angle fit the formula 


sin(2 di) 


2 nit(A t — [/,/ tan/3) 

2 9 

m~ — m~ 

1 1 12 


(29) 


In the sbottom sector on-shell and tree-level masses on the other hand differ by a shift in 
the (1, l)-element, see AAA/ in Ref. [7]. Moreover we employ the scheme which works with a 
dependent bottom-quark mass mb, whereas Ab is defined to be on-shell, see e.g. Refs. [83H86] , 
To allow for maximal flexibility the specification of on-shell squark masses and mixing angles 
is now also possible in case of the MSSM. The Block REN0RMSB0T is not of relevance in such 
input files, since mj, is chosen as dependent parameter, whereas the squark mixing angle Ob 
and the soft-breaking parameter Af } are understood as renormalized on-shell. 

Two options for the pseudoscalar Higgs mixing matrix are accepted as input by SusHi, 
namely the full Higgs mixing matrix, which corresponds to the multiplication 1Z P '1Z G in the 
above notation, but instead also the rotation matrix TZ P can be used as input. Following 
SLHA2 [70] the full matrix (TZ p TZ (: ) l j is provided in Block NMAMIX and asks for entries ij 
with i E {2,3} and j E {1,2,3}. The matrix TZ P can be specified in Block NMAMIXR, which 
only asks for entries ij with {i, j} E {2, 3}. We point out that in contrast to other codes the 
Goldstone boson remains the first mass eigenstate, such that Block NMAMIXR does not ask 
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for entries with i = 1 or j = 1. The elements of the CP-even Higgs boson mixing matrix 
are specified in Block NMHMIX [70]. The Higgs masses need to be given in Block MASS using 
entries 25, 35 and 45 for the CP-even Higgs bosons and 36, 46 for the CP-odd Higgs bosons. 

The block Block EXTPAR still contains the gluino mass as well as the soft-breaking pa¬ 
rameters for the third generation squark sector. Entry 23 for the /r parameter is however 
replaced by entry 65, where the effective value of /r needs to be specified. Moreover entry 61 
asks for the choice of A. SusHi extracts the VEV v s from /r and A. Since the Higgs sectors 
including their mixing are provided, there is no need to provide the parameters k, A k , A\ (or 
mjj±) in the SusHi input, since they do not enter the couplings relevant for Higgs production. 
The Block SUSHI entry 2 specifies the Higgs boson, for which cross sections are requested. 
The CP-even Higgs bosons are numbered 11,12 and 13, the CP-odd Higgs bosons 21 and 
22. Similarly the options 11,12 and 21 also work in the 2-Higgs-Doublet Model (2HDM) and 
the MSSM and 11 and 21 in the SM. A CP-odd Higgs boson 21 in the SM is obtained from 
the 2HDM case with tan/3 = 1. We note that SusHi is still compatible with input files with 
0 (light Higgs), 1 (pseudoscalar) and 2 (heavy Higgs) as options for entry 2. Output files 
however stick to the new convention. 

For the time being we emphasize that SusHi is not strictly suitable for very low values of 
Higgs masses < 20 GeV, where quark threshold effects start to become relevant and also 
electroweak corrections are not implemented. This statement mostly applies to studies of a 
very light CP-odd Higgs boson, which is poorly constrained by LEP experiments in contrast 
to a light CP-even Higgs boson m- 

4 Phenomenological study 

In this section we elaborate on the phenomenological consequences of the additional SU(2)^ 
singlet in the NMSSM with respect to the MSSM for neutral Higgs production. Neglecting the 
squark induced contributions to gluon fusion, the only consequence of the additional singlet 
component is another admixture of the three CP-even/two CP-odd Higgs bosons. However, no 
generically new contributions to Higgs boson production arise. This differs when taking into 
account squark induced contributions to gluon fusion due to the additional singlet to squark 
couplings. In particular for the CP-odd Higgs bosons squark contributions are only induced 
at the two-loop level due to the non-diagonal structure of the CP-odd Higgs bosons to squark 
couplings. Subsequently we work with two scenarios, start with their definition, present the 
Higgs boson masses and admixtures and then discuss the behavior of cross sections, including 
the squark and electroweak corrections to the gluon fusion cross section. Our studies are 
performed for a proton-proton collider with a centre-of-mass (cms) energy of y/s = 13TeV, as 
planned for the second run of the LHC. Lastly we add a short discussion of renormalization 
and factorization scale uncertainties as well as PDF+a s uncertainties for one of the two 
scenarios. 

4.1 Scenarios Si and S 2 

To present the most relevant features of the NMSSM for what concerns neutral Higgs produc¬ 
tion we pick two scenarios. The first scenario Si is in the vicinity of the natural NMSSM [5] 
with a rather large value of A = 0.62. Other input parameters are M\ = 150 GeV, M 2 = 
340GeV, M 3 = 1.5TeV, tan/3 = 2, A K = —20GeV and \x = 200GeV. A\ is determined from 
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the charged Higgs mass m H ± = 400 GeV. The size of A ensures a large mixing of the singlet 
component with the H { ( \ and H® doublets. All soft-breaking masses are set to 1.5 TeV except 
for the soft-breaking masses of the third generation squark sector, which are fixed to 750 GeV. 
The soft-breaking couplings are set to A = 1.8 TeV. The on-shell stop masses are then given 
by = 544.7 GeV and = 941.2 GeV, whereas the sbottom masses are = 749.4 GeV 
and riij n = 757.4 GeV. We vary k between 0.15 and 0.80 and thus vary the mass of the singlet- 
like Higgs component in particular in the CP-even Higgs sector. We note that for illustrative 
reasons the perturbativity limit approximately given by \/A 2 + k 2 < 0.7 is not always ful¬ 
filled in our study. We work out the characteristics for the singlet-like component in the 
following discussion. The relevant input for SusHi is obtained with NMSSMCALC 1.03, which 
incorporates the leading two-loop corrections 0(ot s ott ) to the Higgs boson masses calculated 
in the gaugeless limit with vanishing external momentum [77]. We request NMSSMCALC to 
work with an on-shell renormalized stop sector and add local modifications to the NMSSMCALC 
input routines to read in on-shell parameters rather than DR renormalized parameter^. 
These modifications guarantee identical on-shell stop masses in NMSSMCALC and SusHi. The 
renormalization of the sbottom sector on the other hand is performed SusHi-internally. 

We also choose a second scenario S 2 , in which we vary A to decouple the singlet-like Higgs 
from the Higgs doublets. The detailed choice of parameters is M\ = 150 GeV, M 2 = 300 GeV, 
M 3 = 1.5TeV, tan/3 = 10, A = —2.0TeV, k = 0.2, A n = —30GeV, /x = 130GeV and 
mjj± = 350 GeV. In this scenario we set the soft-breaking masses to 1.0 TeV. The on-shell stop 
and sbottom masses are given by = 824.1 GeV, m^ 2 = 1173.4 GeV, mg = 998.0 GeV and 
mg o = 1008.4 GeV. We vary A between 0.04 and 0.25. For small values of A H\ corresponds 
to the SM-like Higgs boson with mass ~ 121 GeV. The lower bound at A = 0.04 is to 
avoid tiny cross sections for a heavy singlet-like Higgs boson and to keep its mass below the 
SUSY masses thresholds to justify the NLO SQCD expansion employed for the gluon fusion 
cross section calculation. 

Both our scenarios come along with rather light third generation squark masses at the 
low TeV scale. Contrary to the Higgs mass calculations the squark contributions completely 
decouple from Higgs production for heavy SUSY spectra. Our scenarios are chosen to flash the 
phenomenology of an additional singlet-like Higgs boson and thus do not always include a SM- 
like Higgs boson with mass ~ 125 GeV and are partially under tension from LEP searches m 
(for low CP-even Higgs masses below 110 GeV) or LHC searches [8814102] . 

We add for both scenarios the relevant SM input, which includes the MS renormalized 
bottom-quark mass mb(rrib ) = 4.20 GeV, which is translated into a bottom-quark pole mass 
of mb = 4.92 GeV. In SusHi we choose the renormalization scheme, where the bottom-quark 
pole mass enters all occurrences of heavy bottom-quark masses in the loops and the bottom- 
quark Yukawa coupling for the gluon fusion cross section. Bottom-quark annihilation is based 
on the running MS renormalized bottom-quark Yukawa coupling. As pointed out in Ref. [46] 
the gluon densities are hardly dependent on the bottom-quark pole mass fit value of the 
PDF fitting groups, emphasizing that there is no need to adjust the bottom-quark pole mass 
to the PDF fit value for the calculation of the gluon fusion cross section. The top-quark 
pole mass equals m* = 173.3 GeV. The strong coupling constant a s (mz) is set to 0.1172 for 
the calculation of running masses, and is obtained from the corresponding PDF set for the 
cross section calculation. We choose MSTW2008 [ 103 j at the appropriate order in perturbation 


2 We thank Kathrin Walz for instructions on how to modify the NMSSMCALC input routines. 
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theory. Our central scale choices for gluon fusion are m ^/2 for both renormalization and 
factorization scale, and /i° F respectively, and /i° R = and = m^/i for bottom-quark 
annihilation. 

4.2 Higgs boson masses and singlet admixtures 





(b) 



(d) 


Figure 1: (a) Singlet character and (b) masses of the three CP-even Higgs bosons Hi (black), 
H 2 (red, dashed), H% (blue, dotdashed) as a function of At for scenario Si obtained from 
NMSSMCALC 1.03; (c) Singlet character and (d) masses of the two CP-odd Higgs bosons A\ 
(black), A 2 (red, dashed) as a function of At for scenario Si obtained from NMSSMCALC 1.03. 


Subsequently we start with a discussion of the singlet admixture and the masses of 
the three CP-even and the two CP-odd Higgs bosons, which we obtain through a link to 
NMSSMCALC 1.03 as explained beforehand. For scenario Si the singlet component as a func¬ 
tion of k for the three CP-even Higgs bosons is shown in Fig.|l] (a). Clearly, for low values 
of At the lightest Higgs H\ is mainly singlet-like, whereas with increasing At the dominant 
singlet fraction moves from H\ to H 2 and for large values of At to H%. The sum of all singlet 
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components yields Yli = 1- The masses of the CP-even Higgs bosons can be found in 

Fig.[U(b). With increasing k the mass term of the singlet component in gauge eigenstates is 
increasing proportional to kv s , such that the singlet-like Higgs boson can be identified with 
the Higgs boson linearly increasing in mass. Close to k ~ 0.35 H 2 shows the most dominant 
singlet fraction, which will later be visible in the gluon fusion cross section. Scenario S\ 
includes for k > 0.3 a SM-like Higgs boson Hi with a mass of m,H 1 ~ 125 GeV. For very small 
values of k the decay H 2 —» H\H\ opens and leaves a characteristic signature for the SM-like 
Higgs boson H 2 . Note that a light singlet-like Higgs boson lifts the mass of the SM-like CP- 
even Higgs through singlet-doublet mixing, which for our example equals rn jj 2 ~ 153 GeV 
for k = 0.1. The region of small k and a light CP-even singlet-like Higgs boson H\ is largely 
constrained by the LEP experiments [87] . 

Fig.m (c) and (d) show the behavior of the singlet admixture and the masses for the two 
CP-odd Higgs bosons in scenario S\ as a function of k. We point again to the region in the 
vicinity of k ~ 0.35, where the light CP-odd Higgs boson A\ is a pure singlet-like CP-odd 
Higgs boson contrary to the CP-even Higgs boson # 2 , for which H g and components 
remain. The coupling of A\ to quarks vanishes, but the coupling to squarks is still present 
due to the relatively large value of A = 0.62, which will be apparent when calculating the 
gluon fusion cross section. 

For scenario S 2 Fig. [2] shows correspondingly the singlet character and masses for the CP- 
even and CP-odd Higgs bosons. Due to the fixed value of fx = the singlet-like Higgs 

boson increases in mass (proportional to kv s ) with decreasing A and thus for small A H% as 
well as A 2 clearly decouple from the other Higgs bosons. We will later use this setup to show 
the decoupling behavior of the cross sections. Below A < 0.05 both H 3 and A 2 have a singlet 
character, which exceeds > 0.999. 

4.3 Scenario Si: Inclusive cross sections for yjs = 13 TeV 

In this subsection we investigate the gluon fusion a gg and bottom-quark annihilation a bb 
cross sections for scenario Si for y/s = 13 TeV for a proton-proton collider. The subsequent 
statements are however hardly dependent on the cms energy and thus hold for the 7/8 TeV 
LHC runs as well as for more energetic runs. Fig. [3] shows the cross sections for the three 
CP-even Higgs bosons. Naturally the cross sections are strongly dependent on the Higgs 
mass, which are in turn a function of k. Thus, the cross section for the second CP-even 
Higgs bosons H 2 tends to decrease with increasing n. Crucial is the singlet admixture of the 
Higgs boson under consideration. The larger the singlet component |77/j| 2 , the smaller the 
coupling to quarks becomes and thus the more sensitive is the cross section to squark and 
electroweak contributions. For H 2 we observe a cancellation of quark contributions through 
the admixtures with the SU(2 )l doublets around k ~ 0.35, where in turn due to the generally 
small cross section squark but also electroweak corrections to the gluon fusion cross section 
are of large relevance, see Fig.[3] (c) and (d). For small values of k the decay H 2 —>• H 1 H 1 
opens in addition to the large gluon fusion cross section for the singlet-like CP-even Higgs 
boson Hi. The region is therefore constrained by LEP experiments m- Much smoother is 
the behavior for the bottom-quark annihilation cross section, where the direct coupling to 
bottom-quarks is related to the non-singlet character of the Higgs under consideration. In 
an interval around k ~ 0.35 bottom-quark annihilation even exceeds the gluon fusion cross 
section for H 2 despite the small value of tan (3 = 2. 
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Figure 2: (a) Singlet character and (b) masses of the three CP-even Higgs bosons H\ (black), 
H 2 (red, dashed), if 3 (blue, dotdashed) as a function of A for scenario S 2 obtained from 
NMSSMCALC 1.03; (c) Singlet character and (d) masses of the two CP-odd Higgs bosons A\ 
(black), A 2 (red, dashed) as a function of A for scenario S 2 obtained from NMSSMCALC 1.03. 
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(d) 


Figure 3: (a) Gluon fusion and (b) bottom-quark annihilation in pb at yfs = 13 TeV as well 
as (c) squark and (d) electroweak contributions to gluon fusion for the three CP-even Higgs 
bosons Hi (black), H-j (red, dashed), H$ (blue, dotdashed) as a function of k for scenario Si. 
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We show the effect of squark and electroweak contributions to gluon fusion for the three 
CP-even Higgs bosons in Fig.[6] (c) and (d). a g ^ q in Fig.[6] (c) includes stop- and sbottom- 
quark induced contributions at NLO SQCD on top of the quark induced contributions without 
electroweak contributions and compares to the pure quark induced cross section cr gg without 
electroweak contributions. All cross sections include NLO QCD quark contributions and the 
NNLO QCD top-quark induced contributions in the heavy top-quark effective theory. Fig.[3](d) 
accordingly shows the effect of electroweak contributions induced by light quarks following 
Eq. (11511 in combination with Ref. [36] in comparison to the quark and squark induced cross 
section cr^, CD = a gi \ q . Note that in all our figures a gg corresponds to cq^ CD+EW . As expected 
for H 2 the region with small quark contributions induced by the admixture with the H R 
and H R components is in particular sensitive to squark corrections. For the other Higgs 
bosons the squarks corrections in this scenario are incidentally all of the order of 0(— 10%) 
and mostly independent of k. We note that the squark corrections are mainly induced by 
stop contributions, whereas sbottom-induced contributions only account for a small fraction. 
Interestingly, the squark contributions show an interference-like structure with a maximum 
and minimum around k ~ 0.35, whereas the relative electroweak corrections are always 
positive. This can be understood from a sign change in the real part of the quark induced LO 
and NLO amplitude for H 2 at k ~ 0.35, which is of relevance for the squark contributions, 
whereas the imaginary part, more relevant for the electroweak contributions, does not change 
its sign. The size of the electroweak corrections for H 2 follows from a suppression of the 
couplings of the second lightest Higgs H 2 to quarks in contrast to the couplings to gauge 
bosons. Obtaining a pure singlet-like Higgs boson in the CP-even Higgs sector, which neither 
couples to quarks and gauge bosons, rarely happens due to the mixing between both S R 
and H R as well as S R and H R for large values of A. For the SM-like Higgs boson with a 
mass below the top-quark mass the electroweak corrections by light quarks are typically of 
the order of 0(+ 5%) and cover most of the SM-electroweak correction factor. On the other 
hand, for Higgs masses above the thresholds 2> 2 mw or 2 mz the electroweak corrections 
by light quarks are small. The structure visible for H 2 in Fig. [3] (c) for k < 0.3 is induced by 
the thresholds 2 m\y and 2 mz, which the Higgs mass mn 2 crosses between k = 0.1 and 0.3. 
We leave the distortion of distributions, in particular transverse momentum distributions, for 
such a scenario to future studies. 

Similarly we depict the gluon fusion and bottom-quark annihilation cross sections for 
the two CP-odd Higgs bosons in Fig.|3] (a) and (c). The pure singlet-like Higgs boson A\ 
at k ~ 0.352 is clearly apparent, since both cross sections vanish. The corrections through 
squark contributions as shown in Fig. [3] (b) are very large around k ~ 0.352, since squark 
contributions are not suppressed through Higgs mixing, although they only appear at NLO 
SQCD. Electro-weak corrections induced through light quarks are absent for the CP-odd Higgs 
bosons. We note that in the range k = 0.351 — 0.353, where the LO QCD gluon fusion cross 
section for the light CP-odd Higgs boson A\ are tiny, < 10 -5 pb, the prediction for a g ^ q 
with squark induced NLO SQCD contributions for A\ is unreliable, since SusHi calculates 
NLO QCD contributions through the multiplication of one-loop and two-loop contributions, 
where the latter tend to be significantly larger than the former and can thus even induce 
negative cross sections. However, in these regions the tiny cross sections are not of relevance 
for current searches. The cross section (J g ^ q and the relative correction to the vanishing only 
quark induced cross sections of more than 100% need to be taken with care for the CP-odd 
Higgs bosons. The fact that for a pure singlet-like CP-odd Higgs boson the gluon fusion cross 
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section at NLO SQCD completely vanishes due to the absence of a LO contribution motivates 
to take into account NNLO SQCD stop contributions as it was done for the light CP-even 
Higgs boson in Ref. [46] , A first estimate yields tiny, positive cross sections, but we leave an 
inclusion in SusHi to future work. The CP-even Higgs bosons in contrast have a LO squark 
induced contribution, which leaves ag^ q mostly well-behaved. Only in rare cases, where LO 
squark and quark contributions cancel, similar difficulties can arise. 





Figure 4: (a) Gluon fusion and (c) bottom-quark annihilation in pb at y/s = 13 TeV as well 
as (b) absolute value of the relative squark corrections to gluon fusion for the two CP-odd 
Higgs bosons A\ (black), A 2 (red, dashed) as a function of k for scenario S\. 


4.4 Scenario S 2 '. Inclusive cross sections for y/s = 13 TeV 

In this subsection we examine the inclusive cross sections for scenario S 2 , which includes low 
values of A and thus reflects the decoupling of the singlet-like Higgs. The gluon fusion cross 
section together with squark contributions and electroweak corrections through light quarks 
are shown in Fig. [5] for the three CP-even Higgs bosons. With decreasing A the gluon fusion 
cross section for i7 2 /3 rapidly decreases, which has two reasons: First the increasing mass 
naturally decreases the cross section, but second the singlet-like Higgs boson also decouples 
from the other two Higgs bosons - as apparent in Fig. [21 Thus, the indirect coupling to 
quarks is suppressed. Moreover also the direct couplings to squarks are proportional to A 
and thus decrease in size with decreasing A. The relative correction induced by squark 
contributions remains rather constant, see Fig.[5](b). The small interference structure visible 
around A = 0.22 — 0.23 stems from the interchange of the dominant singlet character between 
H 1 and H 2 . Remark that for a SM Higgs the decrease of the gluon fusion cross section 
due to the increase in mass is between uih = 500 GeV and run = 1200 GeV only a factor 
of ~ 47, whereas we observe a decrease of more than five orders of magnitude, thus mainly 
driven by the decoupling. Electro-weak corrections by light quarks as depicted in Fig. [5] (c) 
are completely absent for a heavy singlet-like Higgs boson H 3 , but show a similar pattern 
at large A for H 2 as for small k in scenario Si. The reason is that the Higgs mass m# 2 in 
both scenarios crosses the 2 mz and 2 mw thresholds for the electroweak corrections by light 
quarks. The bottom-quark annihilation cross section shows a similar decoupling behavior 
and is thus not explicitly shown. 

For the CP-odd Higgs bosons we show the corresponding decoupling limit in Fig.[6l where 
the gluon fusion cross section and the bottom-quark annihilation cross section presented in 
Fig. [6] (a) and (c) respectively decrease dramatically for the singlet-like Higgs bosons A 2 . 
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Figure 5: (a) Gluon fusion in pb at yfs = 13TeV, (b) squark corrections and (c) electroweak 
corrections to gluon fusion for the three CP-even Higgs bosons H\ (black), H 2 (red, dashed), 
H 3 (blue, dotdashed) as a function of A for scenario S?. 
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Figure 6: (a) Gluon fusion and (c) bottom-quark annihilation in pb at yfs = 13TeV as well 
as (b) the relative squark corrections to gluon fusion for the two CP-odd Higgs bosons A\ 
(black), A '2 (red, dashed) as a function of A for scenario S 2 . 
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Again the squark corrections shown in Fig. [6] (b) remain rather constant. Given the fact that 
the total cross section for A 2 vanishes with decreasing A, we can therefore conclude that both 
the indirect couplings to quarks for A 2 vanish, but also the direct couplings to squarks vanish. 
Thus, the singlet-like Higgs boson A 2 decouples in the limit A being small. 

We point out that the singlet-like Higgs bosons H 3 and A 2 both approach SUSY mass 
thresholds with decreasing A. Therefore we employ the lower bound of A = 0.04, since for 
lower values of A and thus larger masses of H 3 and A 2 we cannot guarantee the validity of the 
NLO SQCD contributions implemented in SusHi. Ref. [ 46] therefore assigned an additional 
theoretical uncertainty to the heavy SUSY masses expansion. In the decoupling regime we 
checked that the cross sections for Hi, H 2 and A\ coincide with the MSSM cross sections 
obtained for a mixing angle of a = —0.12347 with an accuracy of ~ 10 , which resembles 

the remaining singlet fraction of H 1 , H 2 and A\. 

4.5 Theory uncertainties 

In this section we shortly focus on theoretical uncertainties in the calculation of neutral 
Higgs boson production cross sections. Ref. [46] identified the most important theoretical 
uncertainties for the MSSM, which mostly apply to our discussion of the NMSSM as well. 
Apart from the well-known renormalization and factorization scale and PDF+a s uncertainties 
for cross sections at a proton-proton collider an additional uncertainty for gluon fusion cross 
section is the choice of a renormalization scheme for the bottom-quark Yukawa coupling, 
which is of particular relevance if the bottom-quark loop dominantly contributes. Secondly, 
the fact that NLO SQCD contributions are taken into account in an expansion of heavy 
SUSY masses induces an uncertainty, which grows for larger Higgs masses approaching SUSY 
particle masses thresholds. Thirdly, also relevant for bottom-quark annihilation are missing 
contributions in the resummation A?,, which induce an uncertainty, in particular in the limit 
A;, —>• —1. All of the above theoretical uncertainties as discussed in Ref. [46] apply to the 
NMSSM in a similar way. In contrast to the MSSM however phenomenological studies of the 
NMSSM focus on lower values of tan (3, where both the uncertainty from the choice of the 
bottom-quark Yukawa coupling and the uncertainty induced from unknown contributions to 
A;, are of less importance. A detailed discussion in particular for the singlet-like CP-even and 
CP-odd Higgs boson is left for future work. 

In the following we stick to the commonly studied renormalization and factorization scale 
uncertainties as well as the PDF+a s uncertainties. We present our results just for scenario S 1 , 
since no generically new features appear in other SUSY scenarios. We start with the scale 
uncertainty, where we follow the prescription employed in Refs. [8][46]. We thus consider 
seven combinations of renormalization and factorization scales defined as set C of pairs 
(/x) with /i R = { 777 , 0 /4, rn^/2, m^} and /x F = {m^/A, m^/2, m^} under the constraint 
1/2 < /x fi //x F < 2 for gluon fusion. For bottom-quark annihilation the set is determined from 
/i R = {m 0 / 2 , m^, 2 rritj)} and \i F = { 777 - 0 / 8 ,777.0/4,771^/2} with the constraint 2 < [i R l\i F < 8. 
The minimal and maximal cross sections are obtained according to 

a~ := min {a(/j, R ,n F )} , cr + := max (cj(/x h ,/ i F )} , (30) 

(aoji mf) e C M (maiMf) € 

which we present relative to the cross sections u(/x F ,^ F ) at the central scales /i° R and /i° F . 
They are fj, R = /j,° F = 7710/2 for gluon fusion and g® = 7770 and /i 0 F = m^/A for bottom-quark 
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annihilation. We therefore define the relative uncertainties 




_ = cr - <r(fi° R , fi° F ) 

M ’ <7(Vr, t4) 


(31) 




(a) * (b) 




(c) * (d) 


Figure 7: Scale uncertainties A+ and for (a,c) gluon fusion and (b,d) bottom-quark 
annihilation for (a,b) the three CP-even Higgs bosons H\ (black), H 2 (red, dashed), H 3 
(blue, dotdashed) and for (c,d) the two CP-odd Higgs bosons A\ (black), A 2 (red, dashed) 
as a function of k for scenario Si. 


The scale uncertainties are shown for both the CP-even and CP-odd Higgs bosons in Fig. [3 
for yfs = 13TeV. In case of gluon fusion the SM-like Higgs comes with a scale uncertainty 
of about 0(±1O%) taking into account NNLO QCD top-quark contributions in the heavy 
top-quark effective theory. The scale uncertainty is naturally strongly dependent on the 
individual contributions to the cross section and increases in particular in regions, where the 
top-quark induced contributions are small or quark contributions to the gluon fusion cross 
section cancel and come along with large squark and/or electroweak corrections. The latter 
effect is very pronounced for H 2 in Fig. [7] (a) and for A\ in Fig. [7] (c). 
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For bottom-quark annihilation as shown in Fig. [7] (b) and (d) the scale uncertainty is 
mainly dependent on the Higgs mass, rather than the specific SUSY scenario. The large 
uncertainty for low Higgs masses reflects the need to move toward the four-flavor scheme 
(4FS) [104, 1105 ] in the description of the process. 
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Figure 8: PDF+a s uncertainties for (a) gluon fusion and (b) bottom-quark annihilation for 
the three CP-even Higgs bosons H\ (black), H 2 (red, dashed), H 3 (blue, dotdashed) as a 
function of k for scenario Si. 


As a last step we discuss PDF+a s uncertainties by applying the practical PDF4LHC 
recommendation [10611107] for the MSTW2008 [ 103] PDF sets in order to emphasize the findings 
of Ref. [46j for the case of the NMSSM. For this purpose we combine the results obtained 
with the 41 PDF sets of MSTW2008(n)nlo68cl with the a s uncertainties obtained by the 
PDF sets, which vary a s within the 68% confidence level interval. Fig. [8] shows the PDF+o: s 
relative uncertainties Ap DF+Q with respect to the standard PDF+a s choice for the three CP- 
even Higgs bosons for gluon fusion and bottom-quark annihilation. The standard PDF+a s 
choice equals the zeroth PDF set of MSTW2008(n)nlo68cl together with the standard values 
a s = 0.120 at NLO and a s = 0.117 at NNLO QCD. Similar to the MSSM the uncertainties 
are mainly dependent on the Higgs mass and only slightly dependent on the specific SUSY 
scenario, even for the singlet-like Higgs boson. A very similar result applies to the CP-odd 
Higgs sector and is thus not explicitly presented. It therefore seems sufficient to take over 
the full relative PDF+a s uncertainties from a CP-even or CP-odd SM Higgs boson with the 
same mass, which is easily adjustable to future updates of the PDF4LHC recommendation. 
Taking into account the combination of the newest MMHT2014 [108] . NNPDF 3.0 [ 109 ] and 
CT10 [ 110 ] PDF sets naturally results in larger PDF+a s uncertainties. However, the simple 
recipe to obtain PDF+a s uncertainties just as a function of the Higgs mass is applicable for 
the combination of the PDF sets provided by the PDF fitting groups as well. 
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5 Conclusions 


We presented accurate predictions for neutral Higgs boson production at proton colliders 
through gluon fusion and bottom-quark annihilation in the CP-conserving NMSSM. For gluon 
fusion we adapt the full NLO QCD and SQCD results from the MSSM to the NMSSM, based 
on an asymptotic expansion in heavy SUSY masses for squark and squark/quark/gluino two- 
loop contributions. Top-quark induced NNLO QCD contributions are added in the heavy 
top-quark effective theory. Electro-weak corrections to gluon fusion mediated through light 
quarks are taken into account and the resummation of sbottom contributions for large values 
of tan (5 can be translated from the MSSM to the NMSSM. The latter procedure also applies 
to bottom-quark annihilation. 

Our discussion comes along with an implementation of the neutral Higgs boson produc¬ 
tion cross section calculation in the code SusHi. The Higgs sector (obtained by an NMSSM 
spectrum generator) needs to be supplied through the SusHi input file. We briefly focused 
on the new features of the additional singlet-like CP-even or CP-odd Higgs boson for what 
concerns neutral Higgs boson production. Due to possible cancellations of quark induced con¬ 
tributions, squark and electroweak corrections to gluon fusion can be of greater relevance than 
known in the MSSM, in particular for not too heavy third generation squark mass spectra. 
For a small singlet-doublet mixing term, which can be achieved by lowering the parameter A, 
the singlet-like CP-even and -odd Higgs boson can both be decoupled from the remaining 
MSSM-like Higgs sector. The renormalization and factorization scale uncertainties reflect the 
individual contributions to neutral Higgs boson production in case of gluon fusion, whereas 
scale uncertainties for bottom-quark annihilation as well as PDF+a s uncertainties for both 
production processes mainly remain a function of the Higgs boson mass. 

We leave a more detailed investigation of theory uncertainties to future work. Moreover 
interesting for future studies is an expansion in a light Higgs boson mass rather than heavy 
SUSY masses for what concerns the inclusion of NLO and NNLO SQCD contributions, in par¬ 
ticular since for pure singlet-like CP-odd Higgs bosons NNLO stop-induced contributions are 
the first non-vanishing contributions to gluon fusion. Similarly a discussion of distributions 
and of the necessity of resummation for transverse momentum distributions is timely for the 
real NMSSM, but left for future work. 
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A Formulas: Higgs-squark-squark couplings in the NMSSM 


In this section we present the squark couplings to the five neutral Higgs bosons cj) of the 
NMSSM as implemented in the code SusHi. As noted before, the singlet component S does 
not couple to quarks, such that the couplings of the Higgs bosons to quarks can be taken over 
from the MSSM by replacing the mixing angle a and thus the projection on H ( J and H® by 
the proper mixing matrix elements 1Z S /1Z P . On the contrary the singlet component couples 
to squarks, for which we present the Feynman rules in the form 
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9q,ij 
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with v = l/\/\/2 Gf = v 2 d + v 2 . The couplings g~ - of squarks with indices {i,j} to CP- 
even Higgs bosons with k = (1,2, 3} or CP-odd Higgs bosons with k = (1, 2} are subsequently 
presented in gauge eigenstates 

gf# = 29?,ij + Kk39?,ij ( 33 ) 

g£ij = T^ki • (c/3 g£}j - spg£Jj) + n% 2 ■ {spg£}j + cpgf%) + , (34) 

where in the CP-odd sector the prerotation with 1Z G involving cp = cos j3 and sp = sin/3 
is performed. They were obtained with the code MaCoR Dm and cross-checked against 
the formulas of Ref. |5] for what concerns the LO stop contributions to gluon fusion. The 
individual contributions g p jj in the CP-even sector with k = {1,2,3} yield: 
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In the CP-odd sector contributions with identical squark indices do not exist, 
remaining ones are given by: 
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All occurrences of the soft-breaking parameters and A 5 were replaced by their relation 
to the squark mixing angles 9 b and 6 t■ Trigonometric functions are abbreviated through 
s x = sin x, c x = cos x and t x = tan x. 
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